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1 . Introduction 



The problem of how to utilize extensive test data, accumulated 
under non - envir onmental conditions, to infer results on system re- 
liability under environmental conditions arises during the early stages 
of development of weapons systems. In particular, the problem of 
estimating the ratio of the mean times to failure under environmental 
and non -envir onmental conditions was discussed in a series of reports 
prepared by the Los Altos office of C-E-I-R, Inc. , ( [ 5] , [ 6] , 

and [7] ) for the Special Projects Office. This ratio, usually labeled 
the K-factor, plays the role of a correction factor to transform re- 
sults obtained under non -envir onmental conditions to produce new 
results that would apply under environmental conditions. Since the 
mean times to failure are unknown quantities under environmental as 
well as non -environmental conditions, their ratio, the K-factor, is 
unknown as well. In the reports quoted above, the authors developed 
estimates of the K-factor for the cases when time to failure is meas- 
ured continuously as well as in cycles. They also developed confidence 
interval estimates for the K-factor when time to failure is measured 
continuously. Apparently exact confidence interval estimates of the 
K-factor have not been developed for "cycle-type items", i. e. , when 
time to failure is measured in cycles. Although various approximate 
and asymptotic interval estimates described in the literature 



1 



(see [2] , [4] , and [6] ) are adaptable to this problem, no exact 

confidence interval has been found up to date. 

In this report, we develop a confidence interval for the K-factor 
when cycle-type items are tested. The sampling procedure applied, 
however, is not sampling to failure (as in [ 5] , [ 6] , and [ 7] ) , but 

is instead sampling to a fixed number of successes. 

The confidence interval obtained is one based on the exact 
(conditional) distribution of the observable random variables. The 
confidence interval itself may not be called exact only by virtue of the 
fact that the observable random variables are discrete (see page 513 
of [ 3] ) , but it should be distinguished from confidence intervals based 
on approximate distributions of the observable random variables. For 
this reason, the procedure described in this report is applicable even 
when only a small number of tests may be performed. 
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2. An Exact Confidence Interval for the K-Factor 

We shall assume that a given number of cycle-type items are avail- 
able for testing under environmental and also under non-environmental 
conditions. We shall not assume to have the same number of items 
available for testing under both conditions. In fact, in usual circum- 
stances the number of items available for testing is much larger under 
non-environmental conditions, for obvious reasons. We further 
presume that items are tested a fixed number of cycles, or until they 
fail to function. This procedure is repeated with a given number of 
items under both environmental and non-environmental conditions 
until a fixed number of successes, i. e. , cycles without failure, are 
achieved. The number of failures under environmental and non- 
environmental conditions are two (different) random variables that 
may be observed in the experiments. We introduce the following 
notation. 



Let 



r^ = total number of successes (i. e. , cycles without failure) of all 



items tested under non-environmental conditions; 



r 



2 



total number of successes of all items tested under environmental 



conditions ; 



X 



total number of failures of all items tested under non- 



environmental conditions when tests run to the r 



th 



1 



success; 
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Y = total number of failures of all items tested under environmental 

. . th 

conditions when tests run to the r^ success. 

As the notation suggests, r and r are numbers fixed usually before 

i ^ 

the start of the experiments, and X and Y are random variables ob- 
served at the completion of the experiments. 

To facilitate building our probabilistic model, we stipulate the 
following conditions and notations. 

1) Failure or success during a given cycle may occur independently 
of other failures or successes during other cycles of the same or 
other items. This may restrict the number of cycles an item may 
be tested. Items whose future function is affected by their failure 
during a cycle cannot be tested after their first failure. If repair 
of an item, however, returns it to its original condition, then such 
an item may be tested a fixed number of cycles regardless of any 
number of failures during the tests. 

2) The probability of failure, , during a cycle under non- 
environmental conditions is the same for all cycles of all items. 

The same is true for the probability of failure, , during a 
cycle under environmental conditions. We have then the following 
definitions: 

q^ = probability of failure during a cycle under non-environmental 
conditions ; 

q = probability of failure during a cycle under environmental 
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conditions ; 



= 1 - q = probability of success during a cycle under 
non-environmental conditions; 

p^ = 1 - q^ = probability of success during a cycle under 

environmental conditions. 

3) Success or failure of an item during a cycle under environmental 
conditions does not affect the success or failure of any item under 
environmental conditions. 

The above conditions imply that X and Y are independent random 
variables with the negative binomial distribution of parameters, 
r^ , q^ , and r^ , > respectively. Therefore, 

p«x = k , = v 1 ) p' 1 C 

, r, k “> 

p<Y = k > ■ C 2 V") p 2 2 o 2 

for k = 0 , 1 , 2 , . . . 

Next, we obtain the conditional distribution of X, given the 
condition X + Y = m . 



P(X = k|X+Y = m) 



P (X = k , X + Y = m) 
P (X + Y = m) 



P (X = k) P (Y = m - k) 
P (X + Y = m) 



The distribution of the sum X + Y is computed routinely: 
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P (X + Y = m) 



m 



L P (X + Y = m | X = x) P (X = x) 
x=0 



m 



£ P (X = x) P(Y = m - x) 
x=0 



m . , r x m - x 

r, + x - K + m - x - 1\ 1 2 



■ £ c 1 

x=0 V 



)C Vl 

A m 



x ) P 1 P 2 q i q 2 



Finally, for k = 0 , 1 , . . . , m 



s r | + k 



P(X=k|X+Y = m) = 






1 '\ /• r? + m - k 
) 'V m - k 



*) 



S m (< > : r l ' r 2> 



( 2 ) 



where 



m 



S ( q I r i« r o ) = £ 



z' r j + x - Is r 2 + m - x - 1 \ x 



m 



1 ’ 2 



x=0 



V 



) K 



m - x 



j 



and 



q = 



It is of special significance that the conditional distribution of X, 
given X + Y = m (represented by formula (2) ) , does not depend on 
and themselves but merely on their ratio q. It is well-known 
that under the above conditions the mean number of cycles to the first 



failure (or between consecutive failures) is 



1 



and 



1 



A ^2 

respectively, under non-environmental and environmental conditions. 

^i 

Therefore, the ratio q = is also the ratio of the mean number 
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of cycles to failure under environmental conditions and the mean number 
of cycles to failure under non-envir onmental conditions, i. e. , q is what 
we called the K-factor in the introduction. Since the (conditional) 
distribution of X, given X + Y = m, depends on q alone, a confi- 
dence interval for q may be based on this distribution. Indeed, as in [3], 
we may always find two real numbers q and q such that 



and 



k 

L 

x=0 



g ( x ; q ; r , 
m I 




m 

2 g m (x; q; V r 2 ) 

x=k 



a 

Z 



( 3 ) 



where g (x; q; r , r ) is the conditional probability mass function 
of X , given X + Y = m (presented in formula (2) ) . Then the interval 
( q , q ) is a 100(1 - ot ) % confidence interval for q . This must be 
accompanied with the usual remark that when the observed value of 
X = k = 0 then q = 0 , and when X = k = m (i. e. , Y = 0 ) then 
q = oo , by definition. 

One-sided confidence intervals are obtained in a similar manner. 

If q and q are determined respectively by 

k 

L g m ( X ’ ^ ; r i ’ r 2 ) = “ (4) 

x=0 

and 

m 

L g ( x ; q ; r , r ) = a , (5) 

, m —ll 
x=k 
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then ( q , 00 ) and (0, q) are respectively lower and upper one - sided 
100(1 - a) % confidence intervals for q. The previous remarks con- 
cerning cases when k = 0 or k = m apply here too. 

Unfortunately, the (conditional) distribution g (x ; q ; r , r ) 

m 1 Z 

is not related to any well-known distribution in an explicit way (such 
as, e. g. , the relationship between the binomial and F-distr ibutions) 
to enable us to obtain percentiles of this (conditional) distribution 
without further computations. For this reason, we developed tables 
for the (conditional) distribution g (x; q; r , r ). Some of these 
are presented in tables I and II in the appendix. 
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3. Numerical Examples 



It is to be expected that whenever non -environmental and environ- 
mental test data must be compared, the former will be more extensive 
than the latter, i. e. , in most typical problems we shall find that 

r >> r . For this reason, the two cases given most attention in our 
1 

investigation were: 



case 1) 


r l = 2r 2 


and r 2 = 


300, 


3 50 . 


, 400, 450, 500 ; 


and 


case Z) 


= 10 r 2 


and r 2 


= 10, 


20, 


30, 40, 50 . 





For these two cases, we developed tables for the percentiles of the 
cumulative (conditional) distribution function 

k 

G m (k; q> r i ’ V = 2 g m (X; qi r i ’ r 2> ' (6) 

x=0 

The choice of m was governed by the fact that the number of failures 

under both environmental and non-envir onmental conditions (i. e. , the 

values of Y and X) are usually quite small. For this reason, we let 

0 ^ m ^ 1Z for case 1) and 0 ^ m ^ 10 for case Z) . Tables I and 

II of the appendix list the percentiles of G ( k ; q ; r , r n ) for 

m 1 Z 

cases 1) and Z) , respectively. Other cases of interest have been 
considered. Tables for the distribution function (6) have been prepared 
(but have not been included in this report) for the following additional 
cases: 
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case 3) r = r = 400, 500, 600, 700, 800, 900, 1000 and 
1 2 

0 £ m ^ 10 ; and 

case 4) r ^ = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 and 

0 £ m ^ 20 . 

Next, we present a few examples of confidence intervals for q . 

1) Suppose the experimental data results in four failures under non- 

environmental conditions when tests have been run until the 300^ 

success (i. e. , a total of 304 cycles) , and in three failures under 

environmental conditions when tests have been run until the 30*^ 

success (i. e. , a total of 33 cycles) . Then, r = 300 , r = 30 , 

1 2 

k = 4, and m = 7 . Therefore, from Table II, we find that 
( . 0325 , . 6969 ) is a 90% confidence interval for q . 

2) Suppose for r = 800, r^ = 400, we observed X = 4 and 

Y = 6 . Then, we have k = 4 and m = 10 . From Table I, 

we find that (.0895, 1.151) is a 90% confidence interval for q . 

3) For r^ = 800, r^ = 400, we observed X = 1 and Y = 5 . 

Then, k = 1 and m = 6 . From Table I, we find that 

( . 0004, 1. 474) is a two-sided 99% confidence inter val for q . 

We also find from the same data that ( 0 , . 5241 ) is a one-sided 
90% upper confidence interval for q . 

4) For r - 1000, r^ = 500, we observed X = 2 and Y = 5 . 
Then, k = 2 and m = 7 . From Table I, we find that 

( . 1029 , . 7409 ) is an 80% confidence interval for q . 
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It may be observed that, in general, the confidence intervals for q 
were shorter when k < m - k as compared to when k > m - k . 
Indeed, failures during non -environmental tests are usually less frequent 
than during environmental tests, which indicates that q^ < q^ or 
q < 1 . This fact, of course, is born out by the preceding examples of 
confidence intervals for q . 
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APPENDIX 



TABLE I 



PERCENTILES OF THE DISTRIBUTION FUNCTION 



G (k; r r ) 
m 12 
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TABLE I 

PERCENTILES OF THE DISTRIBUTION FUNCTION G m (k; 
r. =600, r, =300 
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